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One of the greatest uncertainties affecting the design of multiphase flow technologies for space 
exploration is the spatial distribution of phases that will arise in microgravity or reduced gravity. On 
Earth, buoyancy-driven motion predominates whereas the shearing of the bubble suspension controls its 
behavior in microgravity. We are conducting a series of ground-based experiments and a flight 
experiment spanning the full range of ratios of buoyancy to shear. These include: (1) bubbles rising in a 
quiescent liquid in a vertical channel; (2) weak shear flow induced by slightly inclining the channel; 
(3) moderate shear flow in a terrestrial vertical pipe flow; and (4) shearing of a bubble suspension in a 
cylindrical Couette cell in microgravity. We consider nearly monodisperse suspensions of 1 to 1.8 mm 
diameter bubbles in aqueous electrolyte solutions. The liquid velocity disturbance produced by bubbles in 
this size range can often be described using an inviscid analysis. Electrolytic solutions lead to hydrophilic 
repulsion forces that stabilize the bubble suspension without causing Marangoni stresses. We will discuss 
the mechanisms that control the flow behavior and phase distribution in the ground-based experiments 
and speculate on the factors that may influence the suspension flow and bubble volume fraction 
distribution in the flight experiment. 
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Increasing Ratio of Shear to Buoyancy 
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Monodisperse (Potential-Flow) Bubble Suspension d * 1.4 mm 


Electrolytes induce 
hydrophobic bubble-bubble 
repulsion to prevent coalescence 
without Marangoni stresses 



Dual impedance probe: Bubble velocity and volume fraction profiles 
Hot film probe: Liquid velocity 
Video: Bubble size and aspect ratio 
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0.8 M MgS0 4 Increases Viscosity by About 60% 

However, Potential Flow Theory Still Provides 
Accurate Predictions of Drag Coefficient and Aspect Ratio 
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Averaged equations for bubble suspension 

• Bubble phase : 
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Disperse-phase pressure 



Negative pressure due to hydrodynamic interactions leads to 
instabilities on Earth that are absent in microgravity 
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Detection of Instabilities: 
Vertical Channel Studies 


Liquid Velocity Variance 
Much Larger than Expected 
For Homogeneously 
Distributed Potential-Flow Bubbles 



Visual Evidence of Structure: 
Some Horizontal Clustering 
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Instability in Vertical and Inclined Channel 


Frequency Spectrum of Liquid 
Velocity Shows Most of the 
Energy is at Frequencies 
Larger than U/a 


An Instability That is More 
Apparent to the Naked Eye 
Arises at Higher Volume Fractions 
and Inclination Angles 
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Inclined Channel: 

Bubble volume fraction variation 
drives suspension flow 
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Viscosity associated with the instability-induced Reynolds stress 
Is 1 00 times larger than fluid viscosity and 30 times larger than 
viscosity predicted for a homogeneous suspension 



V-[iUe ff V\J] = -pg</> + Vp f 
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Instability induced bubble pressure or diffusivity is also very large 

«(F b +F l +F d ) = VP -► V-[(U B +U L )^]-V-[OV^] = 0 



36 JU f dR d dtp 
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Apparent bubble viscosity and pressure observed 
in ground-based experiments in an inclined channel 
are greatly enhanced by an instability. 

The instability results from the negative pressure due to 
hydrodynamic interactions which would be absent at Og 


NASA/CP— 2004-2 13205/VOL 1 499 


Volume fraction profile in a vertical pipe flow 
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Deficit of bubbles 
near pipe wall 
due to repeated bubble 
bouncing from wall 
which would be absent 


AtOg 



FIG. 12. Sequence of photographs illustrating one cycle of the bubble 
bouncing motion. The bubble radius is 0.7 mm and the inclination angle is 
83°. 
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Couette Flow of Bubble Suspension at Og 


U ~ u L 

Gravity-induced instability absent 
Repeated bouncing from wall absent 
Potential flow approximation 


more accurate in high Re 
microgravity shear flow 
than on Earth 
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Minimum and maximum volume fractions 



Reynolds Number 
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Lattice-Boltzmann simulations for bubble suspension 
at finite Re 

Bubbles modeled as non-deformable spheres 
with no tangential stress boundary conditions 
(specular reflection of lattice gas) 
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Simulation Parameters: 


R/L = 5 
L/d = 7.67 
Re = 0.117 


volume fraction 



fluid and particle velocity profiles 



fluid velocity ■ particle velocity linear profile 
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Conclusions 


A buoyancy-driven instability (not readily apparent to the naked eye) 
greatly enhances the apparent bubble-phase viscosity and pressure 
in an inclined channel flow 

Buoyancy driven bubble-wall interactions create a deficit of bubbles 
near the wall in vertical pipe flow 

These effects should be absent in Og 

We predict that the bubble volume fraction distribution in microgravity 
Couette flow will result from a competition between bubble-phase 
pressure and centrifugal forces 
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